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ABSTRACT 

The CXOCY J220132. 8-320144 system consists of an edge-on spiral galaxy lensing a background 
quasar into two bright images. Previous efforts to constrain the mass distribution in th e galaxy have 



suggested that at least one additional image must be present (Castander et al. 2006). These extra 
images may be hidden behind the disk which features a prominent dust lane. We present and analyze 
Hubble Space Telescope (HST) observations of the system. We do not detect any extra images, but 
the observations further narrow the observable parameters of the lens system. We explore a range of 
models to describe the mass distribution in the system and find that a variety of acceptable model fits 
exist. All plausible models require 2 magnitudes of dust extinction in order to obscure extra images 
from detection, and some models may require an offset between the center of the galaxy and the center 
of the dark matter halo of 1 kiloparsec. Currently unobserved images will be detectable by future 
James Webb Space Telescope (JWST) observations and will provide strict constraints on the fraction 
of mass in the disk. 

Subject headings: gravitational lensing: strong - dark matter - galaxies: structure - galaxies: halos - 
galaxies: spiral 



1. INTRODUCTION 

Gravitational lensing is a potent probe of the matter 
distribution in the central regions of lensing galaxies. A 
particularly promising application of gravitational lens- 
ing is to the study of the mass distribution in edge-on disk 
galaxies. As the geometry of the luminous disk galaxy 
differs significantly from the dark matter halo in which 
the galaxy is embedded, disentangling the relative con- 
tributions to the total mass of the luminous and the dark 
components is, at least in principle, feasible. 

The decomposition of the disk galaxy from the dark 
matter halo has customarily em ployed a maximum disk 



model (van Albada et al. 1985L In systems where the 



mass-to-light ratio of the disk is unknown, the contribu- 
tion of the disk mass to the rotation curve is taken to 
be as large as permitted by the observed rotation curve, 
and the galaxy dominates the mass near the center of the 
galaxy. The accu racy of the maximum disk mod el is un- 
settled, although |Debattista fc Sellwood| ( 1998 1 suggest 
that disk galaxies with bars are likely to be maximal. 
The characteristic radius at which one typically charac- 
terizes the fraction of matter in disk galaxy is 2.2 disk 
scale lengths (2.2Rd), as this is the radius at which the 
circular speed peaks for a thin exponential disk. In what 
follows, we parameterize the size of the galaxy disk by 
the fraction of mass in the disk galaxy inside a sphere of 
radius 2.2Rd, referred to hereafter as fdisk(2-2Rd). 

Progress on this front has until now been impeded 
by the small number of edge-on spiral galaxy lenses. 
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Four lens sy stems with quasar sources are known: 
Q2237+0305 dHuchra et al.|[1985|, B1600+434 (|Jackson 



et al. 19951 [Jau nscn fc Hjorth||l997[ ), PMN J2004 -1349 
dWinn et al.l_.003l) , and CXOCY J220132. 8-320144 (|Cas 



|tander et al.|200 6). In addition the Sloan WF C Edg e-on 
Late-type Lens Survey (SWELLS, [Treu et al.|[201l| has 
assembled 20 galaxy-galaxy lenses. 

As the lens galaxy in Q2237+0305 has a low redshift, 
zi ens = 0.04, the lensing is most ly sensitive to the bulge 
mass (van de Ven et al. 2010 1. Still, combining with 
stell ar kinematic inform ation, Trott & Webster ( 2002 1 
and Trott et al. (2010) find that the lens galaxy has a 
submaximai disk - i.e., the disk does not dominate the 
mass i n th e inner regions of th e galaxy. |Koopmans et aL| 
( |1998| > and Mailer et al. ( 2000 ) put constraints on a iower- 
hmit for the halo axis rati o m B1600+434, fin ding that 
is greater than q ~ 0.5, and |Maller et al. ( 2000[ ) rules out 
a maximal disk. For PMN J2004-1349, onTyi he bulge- 
to-disk mass is constrained ( Winn et al.||2003 ). 

The SWELLS survey team has analyzed one galaxy- 
galaxy lens, SDSS J2141-0001, combining gas and stel- 
lar kinematics data with t he constraints from the lensed 
arc of the s ource galaxy ( Dutton et al.|[201lj |Barnabe| 
et al. 2012). In one analysis, performed before stellar 
kinematic data became available, the disk of the lensing 
galaxy was found to be submaximai, and fdiski^-^Rd) is 
0.45. In the analysis including stellar kinematic data, the 
disk of the lensing galaxy is maximal and fdisk(2.2Rd) is 
0.72. 

CXOCY J220132. 8-320144 (hereafter CX 2201) was 
discovered as part of the Calan-Yale Deep Extra galac- 
tic Research (CYDER) survey of Chan dra fields dCas- 



tander et al. 2003 



Treister et al. 20051. Ground-base' 



optical imaging and spectroscopy were carried out, con- 
firming the lensed nature of the system and measuring 
z = 0.323 and z = 3.903 for the ga laxy and the quasar 



respectively (Castander et al. 20061. At first examina- 



tion, CX 2201 may not appear to be a prime candidate 
for studying edge-on spiral galaxy lenses, since it has 
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Table 1 

CX 2201 Lens Model Constraints 



Figure 1. Combined 160W, 475W, and 814W HST image of CX 
2201. Image A is located above the disk of the galaxy and image 
B is located below. They are separated by ~0'.'8. The size of the 
displayed image area is 9" X 2'.'7. 

only two observed images of the background quasar. It 
displays some intriguing features, however, that make 
it worthy of inquiry. The lcnsing galaxy has no visible 
bulge, further eliminating one possible impediment to a 
clean estimate of the dark matter to disk mass fraction. 
In addition, both lensed images lie close to the plane of 
the disk and have similar fluxes. If the galaxy disk were 
maximal, the matter distribution should be highly elon- 
gated and it might be expected that CX 2201 would be 
a four-image system. Previous efforts to constrain the 
mass distribution in the galaxy have suggest ed that at 
least one a dditional image must be present (Castander 
2006|). These extra images may be hidden behind 
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the disk which features a prominent dust lane. 

In this paper, we present Hubble Space Telescope 
(HST) observations of the system. These were carried 
out at longer wavelengths than the original detection ob- 
servations in an effort to observe extra images that may 
have been extincted by dust in the disk at shorter wave- 
lengths. In addition, we improve upon the astrometry 
and photometry of the observed images and use the im- 
proved constraints to narrow a range of lensing models 
to describe the mass distribution in the system. 

2. HST OBSERVATIONS 

Observations using the HST Advanced Camera for Sur- 
veys (ACS) and the Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) were carried out on 
2006 May 11 and 2006 April 13, respectively (GO:10518). 
ACS observations were taken in the F814W and F475W 
filters; NICMOS observations were taken in the F160W 
filter. 

ACS imaging consisted of 4 dithered exposures taken 
in ACCUM mode in each filter; total exposure times 
were 2228 seconds and 2252 seconds in the F814W and 
F475W filters, respectively. The exposures were reduced 
and calibrated using the CALACS calibration pipeline 
(including the PyDrizzle algorithm). NICMOS imaging 
consisted of 4 dithered exposures taken in MULTIAC- 
CUM mode; total exposure time was 2688 seconds. The 
exposures were reduced and calibrated using the CAL- 
NICA/CALNICB calibration pipeline. Quadrant bias 
was corrected using the pedsub routine in the STSDAS 
IRAF package. 

In each image, the two quasar images straddle the 
galaxy, with the line connecting the two images inter- 
secting the galaxy disk at a point offset from the galaxy 
center. There appears to be no pronounced bulge com- 
ponent in the galaxy. The disk is highly-inclined, but has 
a non-zero projected axis ratio and is thus not perfectly 
edge-on. A dust lane is visible in the ACS images. Fig- 
ure [1] shows a combined 160W, 475W, and 814W HST 
image of CX 2201. 



Parameter 


Image A 


Image B 


Galaxy 


x position (as) a 


-0.225 


-0.210 


0.000 


y position (as) 


0.353 


-0.468 


0.000 


position error (mas) 


0.44 


0.45 


3.57 


relative flux 6 


1.00 


0.72 




axis ratio 






0.12 


disk scale length, R^ (as) c 






1.60 



"Image positions are calculated in the rotated frame where the 
galaxy plane is horizontally aligned and both quasar images are to 
the left of the galaxy center. 
'Relative flux in the F160W filter. 

c Rd increases with decreasing filter wavelength in the three images, 
as expected for a spiral galaxy. R^ is l'.'Ol, 1"57, and 2'.'37 in the 
F160W, F814W, and F475 filters, respectively. 



The imag es were fit using the G ALFIT program ( Peng 
et al. 2002) to extract positions and fluxes of the im- 
ages and the galaxy and other parameters necessary to 
completely characterize the galaxy. The galaxy was mod- 
eled as an exponential disk, while the quasar images were 
modeled as Gaussians (quasi-delta functions). The posi- 
tions, fluxes, and scale radii of the quasar images and the 
galaxy were allowed to vary, as were the galaxy position 
angle and axis ratio. The combined model was convolved 
with the point spread functions (PSFs) of each instru- 
ment, which were calculated using the tinytim program. 
In constructing the PSFs, a frequency power-law of —0.5 
was assumed for the source quasar spectrum. 

The dust lane is quite prominent in the ACS images, 
and is visible even in the NICMOS image. Despite the 
difficulty of fitting with the dust lane obscuring much 
of the galaxy center, the positions of the quasar images 
and the galaxy are consistent in the different filters. The 
results of the best fit are shown in Table [T] where the 
error-weighted average values are given and positions are 
rotated such that the galaxy plane is horizontally aligned. 
The position angle of the galaxy East of North is 140.5°, 
139.3°, and 140.1° in the F160W, F814W, and F475 fil- 
ters, respectively. In every filter observed, the quasar 
images do not straddle the center of the galaxy; instead, 
they are offset from the center of the galaxy and lie left 
of the galaxy center in the rotated frame. This offset is 
a key factor in the mass models used to describe the lens 
system and is discussed further in Section [3j We mea- 
sure an offset of the images from the galaxy center of 
~f//2, which is smaller than the value of ~0'.'3 previously 



measured by Castander et al. ( 2006 ) using ground-based 
data. 

Although GALFIT gives the same galaxy center posi- 
tion in all the filters, the fit is quite uncertain. Indeed, 
measuring the center through other means (such as using 
elliptical galaxy isophotes or masking out the dust lane in 
the fit) gives offsets which vary from ~0 / .'15 to ^0'.'5. In 
addition, the disk scale length, Rd, varied over the differ- 
ent filters, increasing with decreasing wavelength. This 
is consistent with a spiral galaxy with redder old stars in 
the core and bluer star-forming regions in the outer arms. 
However, this variation introduces additional uncertainty 
in defining Rd for models and is discussed further in later 
sections. 

The results of the photometry for each filter are 
shown in Table [5] The flux ratio A /B is higher in the 
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Table 2 

CX 2201 Photometry 



Table 3 

Lens Models 





F160W 


F814W 


F475W 


Image A 


21.67 ± 0.01 


23.46 ±0.01 


25.96 ±0.02 


Image B 


22.02 ±0.02 


23.76 ±0.01 


26.19 ±0.02 


Galaxy 


19.65 ±0.01 


21.19 ±0.01 


22.61 ± 0.01 


Flux Ratio 


1.39 


1.32 


1.23 


(A/B) 









F160W filter than in the F814W and F475W filters; 
that is, image A appears slightly redder than image B. 
Since image A is closer to the disk than image B, this 
suggests differential extinction by a galaxy dust profile 
which decreases with distance from the plane of the 
galaxy. However, other possible explanations for this 
reddening include microlensing and intrinsic variability 
on time scales greater than the time delay between the 



Castander et al. (2006) previously measured a 



images, 
flux ratio A/B of ~1.1. 



3. GRAVITATIONAL LENS MODELING 

We now fit several mass models to the observations. 
The two observed quasar images provide 4 position con- 
straints and 1 flux ratio constraint, a total of 5. For 
models that have additional unobserved lensed images, 
the non-detection of such images gives us additional con- 
straints. 

For every lens model considered, all have common in- 
puts. Flux errors are artificially increased to 20% to allow 
for microlensing or intrinsic variability; likewise, position 
errors are increased to 5 milliarcseconds (mas) to account 
for substructure effects. For every model we test, 2 of the 
free parameters must be used to fit the position of the 
source, x s and y s . A Gaussian prior on the position of the 
center of the lensing galaxy is included with the standard 
deviation given by the measured error in Table [l] 

Sampling the lens model parameter space is performed 
using emcee, the Pytho n ensemble sampling toolkit for 
affine-invariant MCMC ( Foreman-Mackey et al. 2012), in 
order to sample the parameter space of the mass models. 

3.1. Two-Image Lens Models 

Under the assumption that the two observed lensed 
images are the only images created by the lens, the small 
number of observable constraints can overconstrain only 
simple lens models with 4 or fewer free parameters. We 
choose for a simple lens model a singular isothermal el- 
lipsoid (SIE). A SIE is a logical model to test, as many 
strong lens systems with early-type lensing galaxies have 
been successfully modeled by mass profile s with the same 
slope as the SIE ( Koopmans et al. 2006 ) . 



A SIE has a dimensionless projected surface density 

b 

2~C 



«(0 



(1) 



where £ 2 = x 2 + y 2 /q 2 , where q is the axis ratio of the 
mass distribution and in the case where the halo is spher- 
ical and q = 1, b is the Einstein radius and related to 
the 1-d velocity dispersion, a, by b = 4tt(o-/c) 2 Di s /D os , 
where D\ s and D os are the angular diameter distances 
from the lens to the source and from the observer to 
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Parameter 


Description 


SIE 


b 


halo Einstein radius 




1 


halo axis ratio 




X s 


source x position 




Vs 


source y position 


NIE + Chameleon 


b h 


halo Einstein radius 






halo axis ratio 




Sh 


halo core radius 




b g 


galaxy Einstein radius 






source x position 




Vs 


source y position 



the source, respectively. Solutions for the de flection and 
magnification from this model are given by (Kassiola & 
Kovner|1993 Kormann et al.|1994[ ). The model and free 
parameters are summarized in Table [3j 

Models of the two-image system with a SIE have a best- 
fit x 2 ~ 1000. A summary of the results of our modeling 
runs is presented in Table [4] The poor outcomes are due 
to the asymmetry of the observed images. Both images 
are located to the left of the galaxy center (in the rotated 
frame). In addition, they have similar fluxes. The SIE 
model creates solutions where one image is located to 
the left of the galaxy center and the countcrimage (on 
the opposite side of the disk) is found to the right of the 
galaxy center and with a much smaller relative flux. 

A simple solution to this problem would be to shift 
the galaxy center to be collinear with the observed im- 
ages. Such a shift might be justified. The lensing galaxy 
in CX 2201 could be analogous to NGC 4631, a nearby 
edge-on spiral galaxy. NGC 4631 is lopsided, the disk 
of the galaxy extending further from one side of the 
galaxy center than from other other side. In CX2201, 
the dust lane in the galaxy obscures the center of the 
galaxy, and the best-fit image positions vary from 160W 
to 475W and 814W. It is possible that the true center 
of CX 2201 lies offset from the midpoint of the galaxy 
disk. However, NGC 4631 is interacting with another 



edge-on spiral neighbor, NGC 4656 ( Robertsfl968 Web 
achew 1969 Weliachew et al. 1978). in addition, every 



single filter n our data of (JX 2201 shows that the images 
0.2" away from th e galaxy center, and Ma gellan 



observations reported in Castander et al. (20061 show 



the images as ~ 0.3" away from the galaxy center. On 
the whole, it's unlikely that dust lane obscuration is suf- 
ficient to motivate a galaxy center located between the 
images. Adopting such a model for the sake of argument, 
however, achieves good model fits, \ 2 — 1-5. 

3.2. Models with 3 & 5 Images 

A more realistic mass model would include both a com- 
ponent for the dark matter halo and a component to de- 
scribe the disk galaxy. We describe the dark matter by 
a non-singular isothermal ellipsoid (NIE) , which is a SIE 
modified by the inclusion of a core radius, such that the 
dimensionless projected surface density is 

= 2 ( S 2 + £2)1/2 > ( 2 ) 

where s is the core radius. The dark matter halo, then, 
has 3 free parameters associated with it: bh, qh, and Sh- 
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Table 4 

Summary of Lens Modeling Results 



# of Model Images 


Model 


Best-Fit x 2 


Description 


2 


SIE 


1790 


Model puts one image on the other side of the halo center 
instead of having both images on one side. 


2 


SIE - galaxy center offset 


1.5 


Galaxy center offset is not supported by observations. 


3 


NIE + Chameleon - no dust 
extinction 


29 


X 2 for observed images: X^im = 3 - X 2 f° r extra image: 
Xcxtra = ^6. Model parameters are not well constrained. 
Observations of the extra image and a tighter constraint 
on the disk scale length are required to distinguish mass 
models. 


3 


NIE -(- Chameleon — 1 mag 
dust extinction 


12 


x| im = °- 8 - Xcxtra = 11. 


3 


NIE + Chameleon - 2 mag 
dust extinction 


0.5 


X| im = 0.4. X 2 xtra = 0- 


5 


NIE + Chameleon - 2 mag 
dust extinction 


7 


x|im = ^ ■ Xcxtra = 0- Unacceptable fits for observed image 
positions. 


5 


NIE + Chameleon - 2 mag 
dust extinction - halo center 
offset 


0.3 


Xlim = °- 3 - Xcxtra = 0. 


5 


NIE + Chameleon + 1O 9 M 
substructure - 2 mag dust ex- 
tinction 


0.5 


Substructure mass is large; smaller masses do not improve 
fit over no substructure. 











A hard prior is included in the modeling such that b\ x > 0, 
q h = [0, 1] and s h > 0. 

The galaxy disk is described by an exponential profile. 
We approximate the exponen tial profile by the di fference 
of two NIEs, as described by Mailer et al. (2000), where 
the core radii of NIEs are given by s ffi i = {l/3)Rd and 
s 9i 2 = (7/3)i?d- This is referred to as a 'Chameleon' 
profile. Since the axis ratio and the scale length of the 
disk is measured observationally, only 1 free parameter is 
associated with the galaxy: b g . A hard prior is included 
such that b g > 0. 

With the addition of x s and y s to describe the source 
position, and our disk+halo model requires a total of 
6 free parameters. The model and free parameters are 
summarized in Table [3j With 6 free parameters, more 
than 2 images are needed to constrain the model. It 
might be thought that every additional image should 
provide an additional 3 observable constraints. This is 
incorrect for the constraints on non-detections of images. 
There are no real constraints on the location of such im- 
ages. However, there are constraints on the brightness of 
extra images. 

In order to determine how bright an extra image needs 
to be to be detected, we model an extra image by smooth- 
ing the flux in the image over the measured point spread 
function (PSF) and locate it between the two observed 
images, where extra images are expected. Model images 
with fluxes greater than l/10th the brightness of image 
A are clearly detectable in the observations. Figure [2] 
shows that subtracting a model image of that brightness 
is observable as a " divot" in the galaxy disk. We set the 
flux ratio of image A to an extra image of 10 (A/ [extra 
image] = 10) as the threshold of detectability and con- 
sider this a 2.5cr detection. Images of this brightness 
add an additional xcxtra = 6-25 to the % 2 of the model. 
Images that are N times brighter than the cutoff add 




Figure 2. The top panels show the galaxy with images sub- 
tracted, while the bottom panels show the resulting brightness in 
a slice along the plane of the disk. Left: The two bright modeled 
images are subtracted from 160W observation. Center: An image 
modeled with flux 1/20 the brightness of image A is subtracted 
from the left-hand image. Right: An image modeled with flux 
l/10th the brightness of image A is subtracted from the left-hand 
image. 



N x 6.25 to the model. 



A model image could be significantly brighter than this 
cut-off and still be unobservable due to extinction by the 
dust lane in the galaxy. By way of illustration, consider 
the dust extinction towards the center of the Milky Way. 
We observe CX 2201 at a redshift of z = 0.323 at 1.6/xm. 
This corresponds to J band (1.2/xm) observations in the 
rest frame. The dust ex tinction toward the center of the 
Milky Way is A v = 31 dScoville et ; aT| |2003[) and Aj = 8 
( Rieke fc Lebofsky|1985 1. So significant amounts of dust 
extinction might be expected near the centers of edge- 
on spiral galaxies with dust lanes such as CX 2201. For 
this reason, we test models which add extinction to the 
modeled extra image. In a model with 1 magnitude of 
extinction, extra images with a flux ratio of 4 are at the 
limit of detectability. At 2 magnitudes of extinction, the 
limit is a flux ratio of 1.6, and at 3 magnitudes the limit 
is a flux ratio of 0.63. 

3.3. Models with One Extra Image 
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Table 5 

Best-Fit Parameter Values for Model with One Extra Image and 
2 Magnitudes of Dust Extinction 



Parameter Best-fit 68% error interval 



(unmarginalizcd) 



b h 


(") 


0.91 


[0, 8.7] 


Qh 




0.55 


[0,1] 


Sh 


(") 


1.06 


[0, 27] 


b g 


(") 


1.34 


[0.71,2.78] 


x 3 


(") 


-0.13 


[-0.17,-0.09] 


Vs 


(") 


-0.03 


[-0.07,-0.02] 
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Figure 3. Top: The critical curve (thick line) and caustic (thin 
line) of the best-fit model, along with the positions of the source 
(plus marker) and the images (observed: circles, modeled: crosses). 
Bottom: The marginalized probability for the flux ratio of image 
A to the extra image. The left-hand panels show the results us- 
ing no dust extinction; the right-hand panels the results using 2 
magnitudes of dust extinction. 



Using our NIE + Chameleon model, two different lens 
system configurations could give rise to 3 ima ge solutions: 
a mass model without a naked cusp (e.g., Ibata et al. 
1999 Egami et al. 2000 ) and a naked cusp mass model. 
In a model without a naked cusp, any non-singular mass 
distribution should produce an odd number of lensed im- 
ages. Given this, observed lens system 2-image systems 
are truly 3 image systems, where the central image has 
been demagnified sufficiently to be unobservable. The 
predominance of lens systems with even numbers of im- 
ages has been used to limit the core radius in lensing sys- 
tems, as small cores result in the demagnification of the 
central image (Keeton 2003). In this configuration, the 
source lies between the inner diamond-shaped tangential 
caustic and the outer elliptical radial caustic, and the 
two outer images have different image parities. 

Naked cusps occur in systems with highly elongated 
mass distributions, when the tangential caustic extends 
outside the radial caustic. Generically, a source inside 
this e xtension produces three ro ughly collinear images 



(e.g., Keeton & Kochanek 19981. It is commonly re- 



ported that the central image is of similar brightness to 
the outer images. The only current detection of a nake d 



cusp system is APM 08279+5255 ( |Lewis et aL]|2002[ ), 



where the central image is observed to be 0.'2 times the 
brightness of the brightest outer image, although lens 
models put the brightness flux ratio at ~ 0.75. In the 
case of a naked cusp, the two outer images have the same 
image parity. 

A summary of the modeling results using no dust ex- 
tinction and 1 and 2 magnitudes of dust extinction is 
found in Table 0] Here we see that the no extinction 
case results in a poorly-fit model and ~ 2 magnitudes of 
dust extinction is necessary for the modeled image to be 





Figure 4. The 80%, 68%, and 50% contours for the the two- 
dimensional marginalized probability for a model with 2 magni- 
tudes of dust extinction and one extra image. The top panel shows 
that the amplitude of the halo and the disk, b^ and b g , are nega- 
tively correlated to each other. The bottom panel shows that is 
correlated with the size of the core radius in the halo, s^. The disk 
mass fractions are shown in labeled dashed contours; the contours 
shown are } diak (2.2R d ) = [0.3,0.4,0.5,0.6,0.7,0.8,0.9]. 

faint enough to escape detection. The best-fit model with 
one extra image with no dust extinction and the best-fit 
model with one extra image and 2 magnitudes of dust 
extinction are described in Figure [3j which shows the lo- 
cations of the source, the images, the critical curve (in 
the image plane), and the caustic (in the source plane). 
While the caustics in the two best-fit models are differ- 
ent, both are naked cusp solutions. In fact, in every case 
tested, the only acceptable models of the 2 observed im- 
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Figure 5. The marginalized probability for the flux ratio A/ [extra 
image] for a model with one extra image, 2 magnitudes of dust 
extinction, and a smaller disk scale length, Rd = l'.'Ol. 



Figure 6. The probability distribution of the fraction of mass in 
the disk at 2.2Rd for the fiducial model with one extra image and 
2 magnitudes of dust extinction. 



ages in CX 2201 consist of naked cusp solutions or require 
a lopsided galaxy. 

We might expect that the flux of the modeled extra 
image would be dimmer in the case of no dust extinction 
than in the case with 2 magnitudes of dust extinction. 
In fact, while models with more dust extinction have 
better fits, the probability distribution the flux of the 
extra image, as shown in Figure |3j remains the same. 
Whether there is no dust extinction or a large amount of 
dust extinction, the extra image has a flux ratio A/ [extra 
image] ~ 1.9 — 2.4. We describe the best-fit parameters 
below using only the 2 magnitudes of dust extinction 
results as our fiducial result. 

The best-fit model and 68% unmarginalized errors for 
the model with 2 magnitudes of dust extinction are 
shown in Table [5] We display unmarginalized errors 
as they better represent the range of acceptable values 
when we consider the full set of parameters (and not any 
subsets of parameters). The error intervals for the pa- 
rameters are large and the modeling does not robustly 
constrain the parameter values. For example, the en- 
tire permitted range of values for q, the halo axis ratio, 
falls within the unmarginalized 68% error interval. When 
marginalizing over all the other parameters, the proba- 
bility distribution of q is flat. The poor constraints on 
parameter values are due to degeneracies between pa- 
rameters in the model. These degeneracies are illus- 
trated in Figure [4] Parameter degeneracies are best rep- 
resented by a probability distribution for two selected 
parameters and marginalizing over the remaining param- 
eters. Marginalized distributions are narrower than the 
unmarginalized probability distribution cited in Table [5j 
The two-dimensional marginalized probability show that 
the amplitude of the halo and the disk, bh and b g , are 
anti-correlated with each other. In addition, bh is corre- 
lated with the size of the core radius in the halo, Sh, such 
that larger bh is associated with larger Sh- 

In our models we use the mean measured scale length, 
Rd = 1'.'60, but the measured values vary with wave- 
length between R d (") = [1.01-2.37]. We have also tested 
a smaller value of Rd = l'.'Ol. In general, employing a 
different Rd doesn't appreciably change the best-fit pa- 
rameters (not surprising given the size of the unmarginal- 
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Figure 7. The marginalized probability distribution for 
fdisk{2-2Rd) as a function of the distance of the extra image from 
the galaxy center. The contours show where the 80%, 68%, and 
50% of the probability lie. The set of the contours on the left em- 
ploy the fiducial disk scale length, Rd = 1'.'60, while the contours 
on the right use a smaller disk scale length. Rd = l'.'Ol. 

ized errors), but it does have a significant effect on the 
brightness of the extra image found. As shown in Figure 
[5j smaller disk scale lengths are associated with larger 
flux ratios and, therefore, dimmer extra images. 

We measure the fraction of mass in the galaxy disk en- 
closed within a sphere of radius of 2.2 times the disk scale 
length (fdisk(2-2Rd))- The model degeneracies make it 
difficult to constrain parameters, but the probability dis- 
tribution of the galaxy disk mass fraction is not flat. 
Figure [4] shows that the disk mass fraction is strongly 
correlated the parameter that measures the mass in the 
galaxy disk, b g , and weakly correlated with the Einstein 
radius of the halo, bh, because the mass in the dark mat- 
ter halo within 2.2Rd is dependent upon both bh and 
the core radius size, Sh- The value of fdisk for the best- 
fit x 2 parameter values is fdi S k(2-2Rd) — 0.23, indicat- 
ing a sub-maximal disk mass. In Figure [61 we show the 
probability distribution of fdisk(2-2Rd)- The distribution 
is very wide and only very small fractions (< 10%) are 
ruled out. Without additional observational information 
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to break the parameter degeneracies, a definitive state- 
ment about the size of the galaxy disk - whether it is 
maximal or sub-maximal - cannot be made. 

Figure [7] shows the extent to which information about 
the third image can be used to discriminate between the 
model fits. Extra images which are located further from 
the galaxy center are associated with larger fdi S k(2-2Rd). 
The location of the extra image, however, is also highly 
dependent upon the disk scale length, such that smaller 
disk scale lengths are associated with images located fur- 
ther from galaxy center. Constraining the model with 
information about the extra image will require breaking 
the ambiguity surrounding the disk scale length. Now 
let us consider the case of a smaller disk scale length, 
Rd = l'.'Ol. If the extra image is located at d C xtra = 0733 
away from the galaxy center, then the disk mass fraction 
is /disfe(2.2i? ( j) = [0.6—0.7]" within 50% confidence inter- 
vals. If the extra image is located just 20 miliarcseconds 
(mas) further, d ext ra = 0'.'35, then 50% confidence inter- 
val expands all the way to fdisk(2-2Rd) = [0.6—1.0]". So 
milliarcsecond constraints on the position of the extra 
image could provide significant constraints on the disk 
mass fraction. 

The probability distribution of fdi.ski2.2Rd) for the 
smaller disk scale length is tilted to larger values than 
for the fiducial disk scale length. Fractions smaller than 
fdisk{2-2Rd) = 20% are ruled out, and, for the best-fit 
X 2 parameter values, fdi S k(2-2Rd) = 0.34. 



3.4. Models with Mulitple Extra Images 

The parameter space for five-image lens systems also 
consists of naked cusp systems, although the image con- 
figuration is no different from 5-image lenses with buried 
cusps. In this case, the tangential caustic extends out- 
side of the radial caustic, and the source is located within 
both caustics. In general, the five image solution has un- 
observed extra images which lie along the plane of the 
disk. Figure [8] illustrates this image configuration. In all 
cases tested, at least one of the extra images is as bright 
as the extra image in the 3 image configurations. Thus, 
2 magnitudes of dust extinction are necessary in order 
to obscure all the extra images. Regardless the amount 
of dust extinction the modeled positions of the observed 
images (not the flux from the three extra images) is the 
leading factor in the size of the x 2 > and X2im = 7, a sig- 
nificantly poorer fit than the best-fit 3-image solution. 

Good 5-image solutions are difficult to produce because 
no mass configurations exist where two bright images line 
up perpendicular to the plane of the disk and offset from 
the galaxy center. The x position of the observed images 
are x = -0.225 and x = -0.210. The best-fit model 
positions are x = —0.232 and x = —0.200; the model 
tilts the image positions, bringing one closer to the galaxy 
center and one further away. Despite the poorer fit to the 
model, a five-image solution might be preferable because 
very flattened dark matter halos are ruled out. 

One option to find a better fit is to, once again, offset 
the position of dark matter halo. In our two-component 
lens model, we can move the center of the dark matter 
halo while leaving the disk galaxy center at its observed 
values. This may be more plausible than shifting both 
halo and disk to lie collinearly with the lensed images, 
especially as the size of the offset required is small, ~ 1 



Table 6 

Best-Fit Parameter Values for Model with Three Extra Images, 
Dark Matter Halo Offset, and 2 Magnitudes of Dust Extinction 

Parameter Best-fit 68% error interval 
(unmarginalizcd) 
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Figure 8. Left: the critical curve (thick line) and caustic (thin 
line) of the best-fit model with three extra images, a dark matter 
halo offset, and 2 magnitudes of dust extinction, along with the 
positions of the source (plus marker) and the images (observed: 
circles, modeled: crosses). Right: the corresponding probability 
distribution of the fraction of mass in the dark matter at 2.2R^. 

kpc. 

Offsets between the halo center and the brightes t clus- 
ter galaxy have b een observed in galaxy cluste rs (|Allen 
|1998| phan et al7||2010[) and in galaxy groups (van den 



| Bosch et al.| |2005; Skibba et al. 2011). In a sample of 
galaxies from the H I Nearby Galaxy Survey (THINGS), 
the offset between galaxy and halo center was less than 
one radio beam (~ 10" or 150 - 700 pc) for 13 of 
15 galaxies, w ith the two remaining gal a xies exhibiting 



larger offsets ( Trachternach et al. 2008). Kuhlen et al. 



( 2012 1 uses a dark matter plus hydrodynamics simulation 



of the Milky Way (" Eris" ) to show that the dark matter 
halo may be offset from the galaxy center by 300 — 400 
pc, and this offset lies preferentially near the plane of 
the disk. The predicted offset is smaller than necessary 
in the CX 22 01 system, but it is oriented similarily. In 
addition, the Kuhlen et al. (20121 produced such an off- 
set without requiring an external perturber like a passing 
satellite galaxy or an incomplete accretion event which 
might be visible in observations of CX 2201. They sug- 
gest that the offset is produced via a density-wave-like 
excitation by the stellar bar. 

If we offset the center of the dark matter halo to x Q s — 
—0.210 and y a R = 0.000 (leaving the galaxy center at 
the observed location) and with 2 magnitudes of dust 
extinction, we find a best-fit \ 2 = 0.3. The \ 2 from the 
observed images is as good as that found for the 3-image 
best-fit solution. The left-hand panel of Figure [8] shows 
the positions of the images, the critical lines and caustics 
of this best-fit 5-image solution. We show the best-fit 
parameter values and the 68% unmarginalized errors in 
Table [6] The probability distribution of fdisk(2-2Rd), 
shown in the right panel of Figure [SJ rules out very large 
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fractions (> 80%) and disfavors a maximal disk solution. 
The best-fit parameter value gives fdi S k(2-2Rd) = 0.41. 
While the shapes of the distribution is different than in 
Figure |6j the constraints on the disk fraction is consistent 
with the constraints from the 3-image result. 

Another option for improving the model fits for 5 im- 
age lens systems is to include substructure in the lensing 
galaxy. If we introduce a substructure clump that is lo- 
cated in the galaxy disk, on the other side of the images 
from the galaxy center, we can reduce the \ 2 to 0.5. How- 
ever, despite several attempts, only massive substructure 
clumps of ~ 10 9 Mq were sufficient to improve the model 
fit. This size substructure mass would be ~ 10% of the 
disk galaxy mass and too large to be plausible. 

4. DISCUSSION AND CONCLUSIONS 

Edge-on disk galaxies have the potential to disentan- 
gle the relative contributions to the mass from the disk 
galaxy and the dark matter halo, as the geometry of the 
luminous disk galaxy differs significantly from the dark 
matter halo in which the galaxy is embedded. However, 
to date, only a small number of edge-on spiral galaxy 
lenses have been available and have been studied. 

CX 2201 represents a particularly interesting case. 
Consisting of a background quasar lensed by an edge-on 
spiral galaxy into two bright images, the lensed images 
are of similar brightness and lie to one side of the galaxy 
center. Previous efforts to constrain the mass distribu- 
tion in the galaxy have suggest ed that at least one ad - 
ditional image must be present (Castander et al. 2006). 
These extra images may be hidden behind the disk which 
features a prominent dust lane. The image configuration 
is also suggestive of a naked cusp configuration, of which 
only one system has been observed to date. 

4.1. Summary of Results 

We present Hubble Space Telescope (HST) observa- 
tions of the system, better constraining the observable 
parameters of the lens system but without detecting any 
extra images. We have explored a range of models to 
describe the mass distribution in the system, discover- 
ing that a variety of acceptable model fits exist. Our 
conclusions are summarized as follows: 

1. 2-image lensing solutions requires that the galaxy 
center lie between the two observed images - a con- 
straint not supported by current observational ev- 
idence which puts the galaxy center 0.2" (1 kpc) 
away. 

2. Models which create extra images require at least 
one extra image have a flux ratio (image A to extra 
image) of ~ 2 — 3. This is sufficiently bright to be 
detectable, and 2 magnitudes of dust extinction is 
necessary in order to obscure extra images. 

3. All likely solutions with 3 images require a naked 
cusp. However, the parameter space is not well 
constrained. Given the degeneracies in the model 
parameters, it is difficult to put bounds on the 
fraction of mass in the disk within 2.2 times 
the disk scale length. Very small fractions (< 
10%) are ruled out, and the best-fit model finds 
fdisk(2-2Rd) = 0.23. An observational constraint 



on the unobserved extra image would break many 
of the parameter degeneracies and put strict limits 
on the dark matter halo fraction. 

4. The observed images are poorly fit by a 5-image so- 
lution, unless the dark matter halo center is moved 
to lie between the two images. Such an offset might 
be well-motivated by simulations of the Milky Way 
which show a similar offset. Large disk mass frac- 
tions (> 80%) are ruled out and the best-fit model 
finds fdisk{2-2Rd) = 0.41. Given the poor con- 
straints, the disk mass fractions are consistent with 
the 3-image result. 

4.2. Future Observations with JWST 

In order to put robust constraints on the mass distri- 
bution of the lensing galaxy, better constraints on the 
disk scale length and imaging of the extra image or im- 
ages in CX 2201 are needed. In both cases, observa- 
tions at longer wavelengths are the key to reducing the 
impact of dust extinction. The James Webb Space Tele- 
scope (JWST) is expected to observed the near and mid- 
infrared (0.6 — 30/im) with unprecedented resolution and 
sensitivity. In addition to longer wavelengths, JWST ob- 
servations will provide further improvements as it will 
feature a smaller PSF than current HST observations. 

We estimate the amount of dust extinction that may be 
observed in CX 2201 at JWST wavele ngth s using the fit- 
ting formula given by |Calzettilfll997D and |Calzetti et al.| 
(2000). The longest wavelength observation we use is 
at l.o^m H band (1.2/xm in the rest frame). At this 
wavelength our extra images fall below the threshold of 
detection and at least 2 magnitudes of dust extinction 
are required. Let us consider the case where there are 
exactly 2 magnitudes of dust extinction in the if-band. 
If we observe CX 2201 in i^-band (2.2/im), there are 1.1 
magnitudes of extinction and extra images with flux ra- 
tios (A/ [extra image]) smaller than 3.6 are detectable. If 
there are 3 magnitudes of dust extinction in the ii-band, 
then at if-band there are 1.6 magnitudes of extinction 
and extra images with flux ratios smaller than 2.3 are 
detectable. In both these cases, the extra images in our 
fiducial models would be observable (c.f., Figure [3 1. The 
Calzctti relations for dust extinction encompass the far- 
UV to near-IR spectral range and the longest wavelength 
probed is 2.2/xm. If we extend the relation to observa- 
tions in the L band (3.5/im observed, 2.6/xm rest frame), 
just outside of the Calzetti spectral range, we find that 
there are 0.25 and 0.37 magnitudes of dust extinction 
for 2 and 3 magnitudes in the H band, respectively, and 
flux ratios (A/ [extra image]) smaller than ~ 7.5 are de- 
tectable. In fact, if there were 12 magnitudes of extinc- 
tion in our H band observations (more than the 8 magni- 
tudes estimated for the Milky Way), we could still detect 
the extra images predicted by our models in the L-band. 

At longer wavelengths, we would expect the amount of 
dust extinction to be even lower than these limits. How- 
ever, at those wavelengths, attempts to observe extra im- 
ages will be hampered by dust emission, and observations 
of CX 2201 between near-infrared and mid-infrared wave- 
lengths may be where both effects are minimized (see, 

Thus, future obser 



for example, da Cunha et al. 



2008) 



vations with the James Webb Space Telescope (JWST) 
represent the most promising avenue to observing central 
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images and measuring the mass of the disk in CX 2201. 
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